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Aims. We present evidence for finite magnetic helicity density in the heliosphere and numerical models thereof, and 

relate it to the magnetic field properties of the dynamo in the solar convection zone. 

Methods. We use simulations and solar wind data to compute magnetic helicity either directly from the simulations, or 

indirectly using time series of the skew-symmetric components of the magnetic correlation tensor. 

Results. We find that the solar dynamo produces negative magnetic helicity at small scales and positive at large scales. 

However, in the heliosphere these properties are reversed and the magnetic helicity is now positive at small scales and 

negative at large scales. We explain this by the fact that a negative diffusive magnetic helicity flux corresponds to a 

positive gradient of magnetic helicity, which leads to a change of sign from negative to positive values at some radius 

in the northern hemisphere. 
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1. Introduction 

The magnetic field in the heliosphere is a direct consequence 
of the solar dynamo converting kinetic energy of the convec- 
tion zone into magnetic energy. The magnetic field is cyclic 
with a period of 22 years on average, but has also significant 
fluctuations on top of this. These fluctuation can be large 
enough to suppress the number of sunspots to minimum 
levels for decades, for example during the Maunder mini- 
mum. This minimum has been associated with the little ice 
age in the early 17th century, although the solar activity to 
Earth climate relation remains ill understood. Of particular 
interest for space weather are strong variations caused by 
coronal mass ejections (CMEs). These events are believed 
to be a result of footpoint motions of the magnetic field at 
the solar surface, driving strongly stressed magnetic field 
configurations to a point when they become unstable and 
release the resulting energy in an instant. CMEs can shed 
large clouds of magnetized plasma into interplanetary space 
and can accelerate charged particles to high velocities to- 
wards the Earth. The main driver of these ejections is the 
magnetic field, where the energy of the eruption is stored. 

Large-scale dynamos, for example the one operating in 
the solar convection zone, produces magnetic helicity of op- 
posite sign at large and small scale. Here, magnetic helic- 
ity is the dot product of the magnetic field and the vec- 
tor potential integrated over a certain volume. For a long 
time it was believed that CMEs are disconnected from the 
actual dynamo process, but this view has changed in the 
past 10 years. In the regime of large magnetic Reynolds 
numbers, or high electric conductivity, the magnetic he- 
licity ass ociated with the sma ll-scale field, quenches the 
dynamo (Pouquet et al.|, Il976[ ). This is a concept that is 
now well demonstrated usi n^ periodic box sin iulations of 
helically forced turbulence ([Brandenburgl |2QQ1[ ). However, 
astrophysical large-scale dynamos are inhomogeneous and 



drive magnetic helicity fluxes, whose divergence is relevant 
for alleviating what is now often referred to as catastrophic 
quenching of the large-scale dynamo. An important frac- 
tion of these magnetic helicity fluxes is associated with mo- 
tions through the solar surface and their eventual ejection 
into the interplanetary space (|Blackman fc Brandenburgl 
z003j. Connecting the dynamo with the physics at and 
above the solar surface is therefore an essential piece of 
dynamo physics. 

In this paper we review the state of such models and 
their ability to shed magnetic helicity and to produce ejec- 
tions of the type seen in the Sun. We begin by discussing 
a simpler model in Cartesian geometry and turn then to 
models in spherical wedges. Finally, we compare with ob- 
servations of magnetic helicity in the solar wind and dis- 
cuss our finding in connection with earlier dynamo mod- 
els. Full details of this work hav e bee n pub lished else- 
wher e (.Warnecke fc Brandenbura |2Q1Q| : IWarnecke et al] , 
l2Qlll : iBrandenburg et al] , l2Qllf ), but here we focus on an 
aspect that is common to all these papers, namely the na- 
ture of magnetic twist associated with the ejecta away from 
the Sun. 



2. Plasmoid ejections in Cartesian models 

A straightforward extension of dynamos in Cartesian do- 
mains is to add an extra layer on top of it that mimics 
a nearly force- free solar coro na abo ve it. This was done 
by Warnecke fc Brandenburg (|2Q1Q| ) who used a dynamo 
that was driven by turbulence that in turn was driven by 
a forcing function in the momentum equation. To imitate 
the effects of stratification and rotation that are known to 
produce helicity, they used a forcing function that was it- 
self helical. This leads to large-scale dynamos that are more 
efficient than naturally occurring ones that are driven, for 
example, by rotating convection (see, e.g., iKapvla et al.L 
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Fig. 2. Magnetic field structure in the dynamo exterior. 
Field lines are shown in red and the modulus of the cur- 
rent density is shown in pink with semi-transparent opacity. 
Note the formation of a vertical current sheet a b ove th e ar- 
cade. Adapted from IWarnecke fc BrandenburgI (|2Q1Q[ ). 



Fig. 1. Three-dimensional visualization of the magnetic 
field viewed from above. The vertical magnetic field com- 
ponent is color-coded (light/yellow pointing upward and 
dark/blue pointing downward). Note that the field lines 
form a left-handed spiral over the scale of the domain, as ex- 
pected for turbulence with positive helicity at small scales. 
The X-axis points to the ri^ht while the 7y-axis point s up- 
ward. Adapted from IWarnecke fc BrandenburgI (|2QlQf ). 



I2Q1QI : IWarnecke eraP . I2Q12I : iKapvla et all . l2Q12f ). In those 
papers, the produced kinetic helicity is much weaker. As 
mentioned in the beginning, such dynamos can produce 
magnetic fields whose small-scale contribution has magnetic 
helicity of the same sign as that of the forcing function and 
whose large-scale contribution has magnetic helicity of op- 
posite sign. 

In our Cartesian model, the two horizontal directions {x 
and y) are equivalent, but when the large-scale field satu- 
rates, it must finally choose one of the two possible direc- 
tions. This is a matter of chance but, in the case discussed 
below, the field shows a large-scale variation in the y direc- 
tion. The large-scale field settles into a state with minimal 
horizontal wavenumber, which is here (kx^ky) = (0,/ci), 
where ki = 2it/L and L is the horizontal extent of the do- 
main. Note that for fields with variation along the diagonal 
the wavenumber would be \/2 times larger, so such a state 
is less preferred. 

In Figure [T] we show the surface magnetic field of suc h 
a dynamo of the work of IWarnecke fc Brandenbur g^ ('201Q). 
We show color-coded the vertical (line of sight) magnetic 
field component together with a perspective view of field 
lines in the volume above, which we shall refer to as the 
corona region. In addition to fiuctuations, we can see a 
large-scale pattern of the field with a sinusoidal modula- 
tion in the x direction and no systematic variation in the y 
direction. The field lines in the corona region show a spiral- 
ing pattern corresponding to a left-handed spiral. This is 
because the helicity of the forcing function is in this model 
positive (right handed), so the resulting large-scale field 
must have helicity of the opposite sign. 

The large-scale magnetic field is essentially steady in 
the dynamo region and does not change its overall bipolar 
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Fig. 3. Dependence of (J • S)H/(-B^)H,t on time and 
height. Dark/blue stands for negative and light /yellow for 
positive values. For this run the vertical extent of the do- 
main is — Itt < z < -" 
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structure; see Figure [2] for a perspective view. However the 
magnetic field of the corona region shows a time-dependent 
oscillating structure associated with nearly regularly occur- 
ring ejection. The ejection events can be monitored in terms 
of the current helicity, J • B^ where J = V x B / jiq is the 
current density and jiq the vacuum permeability. To com- 
pensate for the radial decline of J • B^ we have scaled it 
by {B'^)H,t^ which is here understood as a combined aver- 
age over horizontal directions and over time. The result is 
shown in Figure [3l 

It is remarkable that the field does not remain steady in 
the outer parts. This can be seen more clearly in a sequence 
of field line visualizations in Figure HI Here, the magnetic 
field is averaged over the x-direction and we show {Bx)x 
color coded together with magnetic field lines as contours 
of {Ax)x in the ?/z-plane. Light/yellow shades correspond 
to positive values, the dark/blue to negative, similar to 
Figure [3l Note that a concentration in Bx emerges from the 
lower region to the outer one. The magnetic field line sur- 
round the concentration and form a shape similar to plas- 
moid ejections, which are b elieved to be a two-d i mensi onal 
model of producing CMEs ('Orto lani fc Schnackl . [19931). At 
a time of t/r = 2881 turnover times, the concentration 
is split into two parts, where the upper one leaves the do- 
main through the upper boundary, while the lower one stays 
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Fig. 4. Time series of the formation of a plasmoid ejection. Contours of {Ax)x are shown together with a color-scale 
representation of {Bx)x] dark/blue stands for negative and light/yellow for positive values, as in Figure [3l The contours 
of {Ax)x correspond to field lines of {B)x in the yz plane. The dashed horizontal lines show the location of the surface 
at z = 0. For this run the vertical extent of the domain is — l^r < z < |7r. 



in the lower layer. Here, r = (lirms^f)""^ is the turnover 
time based on the forcing wavenumber /cf and Urms is the 
root mean squared velocity averaged in the lower layer. At 
t/r = 3041, the field line have formed an X-point in the 
center of the upper layer. In an X-point, field lines recon- 
nect and release large amounts of energy through Ohmic 
heating. In the Sun these reconnection events are believed 
to trigger an eruptive flare, which can cause a CME. A sim- 
ilar behavior can also be seen in more realistic models in 
spherical geometry, as will be discussed in the next section. 



3. Helicity reversals in spherical models 

While similar ejection events are also seen in spherical ge- 
ometry, a more surprising property is a reversal of the cur - 
rent helicity in the radial direction (| Warnecke et al.L l2Qll[ ) . 
In Figure [5] we show visualizations of J • B in meridional 
planes. The Figure shows two examples of coronal ejections 
which we found to be ejected from the dynamo region into 
the atmosphere. At t/r = 358 we can identify a shape that 
is similar to that of the three-part structure of coronal mass 
ejections, observed on the Sun (Low, 1996). It consists of 
an arch of one sign of current helicity in front of a bulk of 
opposite sign and a cavity in between. As show in Figure O 
the ejection leaves the domain on the radial boundary and a 
new ejection of opposite sign occurs. We also see that in the 
dynamo region the current helicity is negative in the north- 
ern hemisphere and positive in the southern. This seems to 
be basically true also in the immediate proximity above the 



surface, but there is now an increasing tendency for the oc- 
currence of magnetic helicity of opposite sign ahead of the 
ejecta. This seems to be associated with a redistribution of 
twist in the swept-up material. The current helicity is by 
far not always of the same sign, but both signs occur and 
there is only a slight preference of one sign over the other. 
This is seen more clearly in Figure [6l where we show a time 
series of J B versus radial position r/R and time t/r. 

A time series of the normalized current helicity, 
J ■ B/{B'^)t^ evaluated at radius r = 1.7 R and 28° lati- 
tude, is shown in Figure [71 where we also show their run- 
ning means. It is now quite clear that on average the sign of 
current helicity has changed relative to what it was in the 
dynamo region. This is seen explicitly in a time-averaged 
plot of current helicity in the meridional plane (Figure [8]). 

This reversal is significant because similar behavior has 
also been seen in recent measurements of the magnetic 
hehcity spectrum in the solar wind (Brandenburg et al., 
2011), but before showing the evidence for this, let us first 
discuss in more detail a similar diagnostics for the simu- 
lations. In Figure [9] we show the time variation of Bg and 
^0 and the magnetic helicity spectrum obtained from these 
time series. Unlike the solar wind, where a time series can 
be used to mimic a scan in distance space (the so-called 
Taylor hypothesis), this argument fails in the present case, 
because no wind is produced in the present model. There is 
only the pattern speed associated with the CMEs. Whether 
this is enough to motivate the use of the Taylor hypothe- 
sis is rather unclear. Nevertheless, there is a remarkable 
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Fig. 5. Time series of coronal ejections in spherical coordinates. The normalized current helicity, J • B/{B'^)t, is shown 
in a color-scale representation for different times; dark/blue stands for negative and light/yellow for positive values, as 
in Figure [3l The dashed lines show the location of the surface at r"^ = x'^ -{- z'^ = R'^ . 
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Fig. 6. Dependence of the dimensionless ratio J • B/{B'^)t 
on time and radius. The top panel shows a narrow band in 
in the northern hemisphere and the bottom one a nar- 
row band in the southern hemisphere. Dark/blue stands for 
negative and light /yellow for positive values. The dashed 
horizontal lines show the location of the surface at r = i?. 

similarity with similar helicity spectra obtained for the so- 
lar wind; see Figure [TOl In both cases, magnetic helicity 
has been obtained under the assumption of local isotropy 
of the turbulence. This means that one computes the one- 
dimensional magnetic energy El^{kR) and magnetic hehc- 
ity spectra H^{kR) simply as 



Euikn) = I^IVmo, i^M (^i?) = 4lm{BeB;)/kR 



(1) 



where kR is the component of the wave vector in the radial 
direction. Here, a hat denotes Fourier transformation and 



Fig. 7. Dependence of the dimensionless ratio J • B/{B'^)t 
on time at radius r = 1.7 R and 28° latitude. The solid line 
stands for the northern hemisphere and the dotted for the 
southern hemisphere. The red lines represent the cumula- 
tive mean for each hemisphere. 



an aste risk complex conjugatio n. These are the equations 
used by iMatthaeus et al.l (|l982l ) who applied such an anal- 
ysis to data from Voyager 2. Since Voyager 2 flew close to 
the ecliptic, the magnetic helicity is dominated by fluctu- 
ations. This is why [Brandenburg et all (|2Qllf ) applied this 
analysis to Ulysses data, where a net magnetic helicity was 
seen for the first time. An important advantage of Ulysses 
over Voyager 1 and 2 is the high angle with the ecliptic. 
So, only with Ulysses we can measure the magnetic helic- 
ity far away from the ecliptic in both hemispheres of the 
heliosphere. 
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Fig. 8. Current helicity averaged over time. Dark/blue cor- 
responds to negative values, while the light/yellow corre- 
sponds to positive value, as in Figure [6l The dashed line 
show the location of the surface at r^ = x^ + z^ = i?^. 
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Fig. 9. Helicity in the northern outer atmosphere. The val- 
ues are written out at the point, r = 1.5 i?, 90° — ^ = 17°, 
and = 9°. Top panel: Phase relation between the toroidal 
B(fy and poloidal Bq field, plotted over time t/r. Bottom 
panel: Helicity H{k) is plotted over normalized wave num- 
ber kR. The helicity is calculated with the Taylor hypoth- 
esis using the Fourier trans formation of the poloi dal and 
toroidal field. Adapted from I Warnecke et al.l (|2Qllf ). 



What we see both in the simulations and in the solar 
wind is that there is magnetic helicity of opposite sign at 
large and small scales. However, it is exactly the other way 
around than what is found in the corona region; see Table [H 
To understand the reason for this, we need to consider the 
equations for the production of magnetic helicity at large 
and small length scales. As will be argued in Section SI 
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Fig. 10. Magnetic energy and helicity spectra, 2/io^M(^) 
and k\HM{k)\^ respectively, for two separate distance in- 
tervals (first and third panels). Furthermore, both spectra 
are scaled by AttR'^ before averaging within each distance 
interval above and below 2.8 AU, respectively. The relative 
magnetic helicity, kHM{k)/2fioEM{k)^ is plotted separately 
(second and fourth panels) together with its cumulative 
average starting from the low wave nu mber end. The zero 
l ine is shown in dashed. Adapted from iBrandenburg et al.l 
(I2011D . 



the mechanism that sustains negative small-scale helicity in 
the north is turned off in the solar wind, and there is just 



Table 1. Distribution of magnetic helicity at large and 
small scales both in the dynamo region close to or below 
the surface and the corona region, solar wind or the exterior 
of the dynamo. 



Magnetic helicity 



Large scales small scales 



Dynamo region, interior 

Corona region, solar wind, exterior 
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the effect of turbulent magnetic diffusion which contributes 
with opposite sign. By contrast, inside the dynamo region, 
turbulent diffusion is subdominant, because otherwise no 
large-scale magnetic field would be generated. However, in 
the wind we do not expect the dynamo to be exited, so here 
diffusion dominates. 



4. Connection with earlier dynamo models 

The purpose of this section is to make contact with dy- 
namo theory and to understand more quantitatively why 
the magnetic helicity reverses sign with radius. In essence, 
we argue that the profile of magnetic helicity density must 
have a positive radial gradient to maintain a negative dif- 
fusive magnetic helicity flux and that this is the reason for 
the magnetic helicity to change from a negative sign to a 
positive one at some radius in the northern hemisphere. 

We begin by discussing the magnetic helicity equa- 
tion. The magnetic helicity density is h = A • B, where 
S = V X A is the magnetic field expressed in terms of 
the magnetic vector potential A which, in the Weyl gauge, 
satisfies dA/dt = U x B — riji^J . The evolution equation 
ioT A • B is then 



d_ 
dt 



A B = -2r]fioJ B -V F, 



(2) 



where F is the magnetic helicity flux. Next, we define large- 
scale fields as averaged quantities, denoted by an overbar, 
and small-scale fields as the residual, denoted by lower case 
characters, so the magnetic field can be s^lit into two_con- 
tributions via B = B-\-b. Likewise, A = A + a, J = J-\-j^ 
and U = U -\-u. The evol ution e quation for the mean mag- 
netic helicity density h = A • B is given by 



d_ 
dt 



A B = -2r]iiQj B -V F. 



(3) 



To determine tjie magnetic helicity density of the mean 
field, hm = A- B, wejise the averaged induction equation 
in the Weyl gauge, dA/dt = UxB-\-uxb — tj/j^qJ^ so 
that 



d_ 
dt 



A B = 2uxb B - 2r]i^oJ B -V F^ 



(4) 



The magnetic helicity equation for the fluctuating field, 
hf = h — hm = ci ■ b^ takes then the form 



d 

7— a b = -2u X b B 

dt 



2r]fioj-b-V -Ff, 



(5) 



so that the sum of Equations ^ and ([5]) gives Equation (J3j) . 
Here, the total magnetic helicity flux consists of contribu- 
tion from mean and fluctuating fields, denoted by subscripts 
m and f, respectively, i.e., Fm + i^f = F. Note that, even in 
the limit 77 ^ and in the absence of fluxes, magnetic he- 
licity at large and small scales is not conserved individually, 
but there can be an exchange of magnetic helicity between 
scales. 

A note regarding the gauge-dependence is here in order. 
Obviously, Equation (J5j) depends on the gauge choice for a. 
However, if we are in a steady state, and if hf also happens 
to be steady (which is not automatically guaranteed), then 
we have 



and since the right-hand side of this equation is man- 
ifestly gauge- invariant, V • Ff must also be gauge- 
inva riant. T his pro perty was used in earli e r wo rk of 
Mit ra et al.l (|2Q1Q|), iHubbard fc Brandenburd (|2QlQf ) and 
I Warnecke et al.l (|2Qllf ) to determine the scaling of F^ with 
V/if and thus the turbulent diffusion coefficient hZh- In ad- 
dition, if there is sufficient scale separation between large 
and small scales, which is typically the case in the non- 
l inear regime at th e end of the inverse cascade process 
(Brande nburgl |2001[ ) , then h{ can be expressed as a density 
of linkages, which is itself manifest ly gauge-independent 
(|Subramanian fc Brandenburgl l2QQ6f ). This property then 
also applies to the flux Ff . 

The correlation u xb is known to have two contribu- 
tions, one proportional to B with a pseudo-tensor in front 
of it (the a effect, responsible for large-scale fleld gener- 
ation), and one proportional to J with a coefficient rjt in 
front of it that corresponds to turbulent diffusion, i.e., 

u X b = aB — tjiIjlqJ^ (7) 

where we have again assumed isotropy. The reason why 
the mean magnetic helicity density of the small-scale field 
is negat ive in the north is because a > in the north 
(e.g. Kr ause and Radlerl . 1980), producing therefore nega- 
tive magnetic helicity at a rate —2aB'^ < for small-scale 
fields and -^2aB'^ > for large-scale fields, so that their 
sum vanishes. There is also turbulent magnetic diffusion 
which reduces this effect, because r^t > and J B > 
in the north. In the solar wind no new magnetic field is 
generated, so turbulent magnetic diffusion could now dom- 
inate and might t hus explain a rev ersal of magnetic helicity 
density (Branden burg et aI]j2Qll^ . 

Support for a reversal of the sign of magnetic helicity 
was first seen in dynamo simulations with magnetic helicity 
fiux in the exterior. In Figure[TT]we show a representation of 
magnetic helicity density of small-scale fields h f = a • b ver- 
sus z and t for a model similar to that of Brande nburg et aP 
(2009), but where the magnetic helicity flux is caused by a 
wind that is then running into a shocHj where the flux is 
artiflcially suppressed at height z = 3H. This flgure shows 
that there is a clear segregation of negative and positive 
small-scale magnetic helicity in the dynamo regime and the 
exterior, respectively. 

In our description above we have suggested that the 
magnetic helicity production balances the da ■ b/dt term, 
but this cannot be true in the steady state. Instead, 
it m_u_st be the divergence o^the magnetic helicity flux, 
V • i^f. Let us assume that ^f can be approximated by 
a Fickian diffusion law, i.e., Ff = —K h^hf. Simulations 
have suggeste d that tZh /vt is aroun d 0.3 (iMitra et al.U201Ql : 
IHubbard fc Brandenburd . I2OIOI : IWarnecke et al.l . I2OIID . 
Thus, balancing now the source S{z) = —2aB'^ -h 
2r]t/jjoJ • B against the divergence of the flux of magnetic 
helicity at small scales, we have, in a one- dimensional model 
(neglecting the molecular diffusion term, 2r]jj.oj • 6): 



S{z) = -Kh 



dz2 • 



(8) 



V ■Fi = -2uxb-B - 2r]/j.oj • 6, 



(6) 



Taking as an example a source where S{z)/nh = —1 in 
the dynamo interior {z < 0), and S{z) = in z > 0, 

^ Note that the shock at z/H — 3 becomes eventually under- 
resolved and the simulation has to be terminated. This is what 
causes the wiggles in the proximity of the shock. 
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Fig. 11. zt diagrams of hf for an a^ dynamo with a wind 
which stops all of a sudden at z/H = 3. The white hori- 
zontal line marks the location z = H. Light/yellow shades 
indicate positive values and dark/blue shades indicate neg- 
ative values. 



we have a family of solutions of Equation (|8]) that only 
differ in an undetermined integration constant correspond- 
ing to a constant offset in the flux; see the second panel 
of Figure [121 The solutions for which the magnetic helic- 
ity flux, — /^/id/if/dz, is negative in the exterior are those 
for which /if reaches an extremum below the surface. This 
seems to be what happens both in simulations and in the 
solar wind. We can thus conclude that the reason for a sign 
change is not a dominance of turbulent diffusion in the so- 
lar wind, but just the possibility of the magnetic helicity 
density reaching an extremum below the surface (dashed 
red and solid black lines in Figure [12]) , not at the surface 
(dotted blue lines). 



5. Conclusions 

As we have seen in the present paper, magnetic twist (or 
helicity) plays an important role for the solar dyn amo 
(jBrandenburd . I2QQ1I : iBlackman fc Brandenb"ur5j 2QQ2l and 
for p roducing eruptions of the form of CMEs ([L ow, 1996, 
I2QQ1I). The r ecent wor k oflWarnecke fc Br andenburgI ([20101) 
and IWarnecke et alj (|201lf ) tries to combine both aspects 
into one. Although the models are still rather unrealistic 
in many respects, they have already now led to useful in- 
sights into the interplay between dynamo models and so- 
lar wind turbulence. In particular, they have allowed us 
to understand the properties of magnetic helicity fluxes. 
We have confirmed that the hemispheric sign rule of mag- 
netic helicity does not extend unchanged into the inter- 
planetary space, but we have now shown that it must flip 
sig;n somewhere abov e the s olar surface. On the other hand, 
iBothmer fc SchwennI (|1998[ ) found that the magnetic clouds 
follow Hale's polarity and that the sign of the magnetic he- 
licity is the same as in the interior. However, this result is 
not based on rigorous statistics. 

Future work in this direction should include more re- 
alistic modeling of the solar convection z one. Preliminary 
work in this direction is already underway (| Warnecke et al.L 
|2Q12[ ). Furthermore, it will be necessary to allow for the de- 
velopment of a proper solar wind from the dynamo region. 
One of the difficulties here is that, if the critical point is 
assumed to be too close to the solar surface, which would 
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Fig. 12. Sketch showing possible solutions hi{z) (upper 
panel) to Equation (|8]) with S = const = — 1 in z < 
and S = in z > 0. The red (dashed) and black (solid) 
lines show solutions for which the magnetic helicity flux 
{—tvh^h^/dz^ see lower panel) is negative in the exterior. 
This corresponds to the case observed in the Sun. The blue 
(dotted) hue shows the case, where the magnetic helicity 
flux is zero_above the surface and therefore do not reverse 
the sign of hf{z) in the exterior. 



be computationally convenient because it would allow us to 
use a smaller domain, the mass loss rate would be rather 
high and could destroy the dynamo. In future work we will 
be trying to strike an appropriate compromise that will al- 
low us to study the qualitatively new effects emerging from 
this. 
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